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INTRODUCTION 

The investigation of the distribution of energy in the flow of 
a compressible gas that moves without resistance in a straight 
cylindrical pipe and gives off heat indicates that the existence of 
a "thermal" Laval nozzle in the supersonic region is possible 
(reference 1) (that is, B a > 1, the cooling of the flow under 
these conditions results in a rise in the Bairstow (Mach) number). 
This interesting result is actually strongly distorted by the effect 
of the resistance. For the simple case of gas cooling (heat con- 
duction at the wall), the well-known relation between the resistance 
and the heat transfer makes the existence of a "thermal" Laval 
nozzle improbable. If only the heat transfer by contact is taken 
into account, the existence of a "thermal" nozzle is impossible; 
however, if radiation from the products of combustion is also 
considered, the "thermal" nozzle is possible only in a narrow range 
of high temperatures and for large dimensions of the nozzle (pipe 
diameter) . 


1. SOLUTION OF EQUATIONS 


Assumptions are based on the analysis of the simultaneous 
solution of the following equations: the momentum equation (in 

which the friction forces are taken into account by the usual 
"hydraulic" resistance coefficient) 
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dw _ _ 1. dp y w_ 
dx p dx * 2D 


( 1 ) 


the equation of continuity 


wp = constant 


( 2 ) 
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and the heat -transfer equation 




( 3 ) 


where v (m/sec) is the velocity of the flow at a certain section of 
the pipe; p (kg/m^) is the static pressure at the same section; 

T (°K) is the temperature of the gas; p ( (kg) (sec^)/m^) is the 
density; D (m) is the diameter of the pipe; £ is the coefficient 
of resistance; G (kg/sec) is the weight of the gas per second; 

Cp (cal/(kg)(°C)) is the specific heat of the gas; a (cal/(m2)(hr)(°C)) 
is the heat transfer coefficient; Ba = w/c is the Bairstow number 
O&ch number); c (m/sec) is the local velocity of sound; t Q (°K) 

is the temperature of the wall; and dS is an element of the pipe 
area that corresuonds to an element of the length dx. 


The solution of this system for the general case was given by 
the author and M. D. Weisman in 1934 (reference 2) and leads to 
complicated expressions not capable of a clear qualitative analysis. 
In investigating the problem of the formation of a "thermal" Laval 
nozzle, the system was reduced by the author to nondimens ional 
variables and was solved under the assumption t Q = 0 (that is, 

at a wall temperature negligibly small in comparison with the 
stagnation temperature of the flowing gas). 

The conclusions as to the possibility of formation of a 
"thermal" Laval nozzle for tg = 0 will naturally be the more 
favorable in this sense because for certain wall temperatures 
comparable with the gas temperature, the intensity of cooling 
will be less than in the case considered by the author. 

The transformed system of equations is written as follows : 

l a(pe) = _ JL b 

e d£ kn 2d 5 (1) 

( 2 ' ) 


( 3 ’) 


NACA comment; The symbol t^ in equation (3) appears subse- 
quently as tp. Brackets and braces in equations (3), (1*), 

and (3*) do not appear in the original version but are con- 
sidered desirable in the interest of clarity. 


= constant 

e 


til [ 1 +/¥) 4 ) 
9 [ 1+ (¥) f 
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where 8 = t/Tq la the nondimens lonal temperature (Tq is the ini- 
tial gas temperature); it = p/Pq ia nondimens ional pressure 
(p 0 is the initial pressure); M * a/wgc p is the Margulis 
criterion {(g) (m) /sec 2 ) is the gravitational acceleration); 

3 = Ba 2 ; and £ = x/D is the nondimens ional length. 


By elementary transformations , the system (1')# ( 2 * ) , 
(3') is reduced- to the equation 


-f] * & #)■»] 

d £ 1-3 


and 


(4) 


that for constant values of £ and M may he integrated^. By 
simple transformations, corresponding expressions for 8 and it 
can also he written. The relations (4'), (4''), and (4 1 ' 1 ) are 
the result of the integration 2 
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The constant of the values £ and M and the existence of the 
normal relations between them for the case of large values of 
Bairstow (Mach) number has been confirmed several times in both 
Russian and other foreign literature (references 3 and 4). 

p 

Hereinafter \ = (k-l)/2. 
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2. ANALYSIS OF SOLUTION 

(a) Radiation neglected. The hydrodynamic theory of heat 
exchange gives the relation "between £ and M 

£ = 8M (5) 

By substituting equation (5) in equation (4), the following expression 

'Z 

is obtained : 

n - £■(?-}) < 3 * * 6) 

A study of equation (6) shows the existence of three regions 
of variation of (3 : 

1* 3 > 1. In this casej p drops along the pipe when 

approaching p = 1. 

2. l/k < p < 1. Here p increases when approaching p = 1. 

3. p < l/k. In this region p drops and approaches the normal 
flow of an incompressible gas in the presence of cooling. The 
behavior of 3,0, and n for particular cases is shown in 
figures 1, 2, and 3. 

(b) Radiation considered. In the case of the presence 

of triatomic products of combustion (water vapor and carbon dioxide) 
in the flowing gas, a certain quantity of heat will be given to the 
walls by radiation. This condition leads to an increase in a and 
a breakdown of relation (5). Determination of the increase in a 
that is required so that a "thermal" nozzle may exist in region 1 
is made possible by equation (6). It is thus necessary that 



3 

The application of relation (6) facilitates finding the heat stress 

of the endothermal reaction and makes the existence of a thermal 

nozzle fundamentally possible. 
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that is 


M> -JL (7) 

r*nrr ' J 


An increase in the value of a by 17 percent corresponds to equa- 
tion (7) as compared with the value given ty the hydrodynamic 
theory of heat Interchange. 

The computations made for a nozzle with a 200-millimeter 
diameter and hy using the air products of gasoline combustion with 
an excess coefficient equal to 1 indicate that for Pq = 1 at 
atmosphere, the coefficient of heat transfer increases by 
10 percent because of the radiation. In working with a larger 
air-fuel ratio or oxygen-fuel ratio, compositions for which a 
increases by more than 17 percent are possible. This increase 
in a will however occur over a small part of the pipe. After 
a certain lowering in the temperature and a corresponding decrease 
in the radiation, a drop in Ba along the pipe begins (as in the 
usual case) . 1 


3. CONCLUSION 

1. The hydrodynamic theory of heat exchange applied to the 
investigation of the possibility of the formation of a "thermal" 
Laval nozzle on cooling the gas by heat conduction at the wall 
leads to a negative result; that is, the formation of such a 
nozzle is impossible. 

2. When radiation is considered for the case of the flow of 
gasoline combustion products in an air- rich mixture, a certain part 
of the pipe in the region of high temperatures may work as a thermal 
nozzle. 

After a certain lowering of the gas temperature, however, the 
pipe will operate normally (with a drop in the Bairstow (Mach) 
number along the pipe) . The effect of an increase in Mach number 
in this case holds only for: (a) relatively large diameters of the 

pipe and (b) products of combustion obtained in an air -rich mixture. 
The practical application of the thermal nozzle even under these 
conditions is in the author's opinion impossible. 

Translated by S. Reiss, 

National Advisory Committee for Aeronautics . 
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Figure 2. - Dependence of — , n f and 0 on \ for £ =0.016 and p 0 = y . 
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